Draft version January 20, 2013 

Preprint typeset using 1^1^^ style emulatcapj v. 03/07/07 



THE SPATIAL AND VELOCITY DISTRIBUTIONS OF HYPERVELOCITY STARS 

FUPENG ZhANG"*^, YOUJUN Lu^ , AND QlNGJUAN Yu^ 
^ National Astronomical Observatories, Chinese Academy of Sciences, Beijing, 100012, China; luyj@nao.cas.cn 
^ Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing, 100871, China; yuqj@kiaa.pku.edu.cn 

Draft version January 20, 2013 

ABSTRACT 

Hypcrvclocity stars (HVSs) found in the Galactic hafo arc believed to be the dynamical products 
of interactions between (binary) stars and the massive black hole(s) (MBH) in the Galactic center 
(GC). In this paper, we investigate how the spatial and velocity distributions of HVSs are connected 
with their originations, ejecting mechanisms, and the dynamical environment in the GC. It has been 
shown that the detected HVSs are spatially consistent with being located on two thin disks (Lu ct 
al.), one of which has the same orientation as the clockwise rotating stellar disk in the GC. Here, we 
perform a large number of three-body experiments of the interactions between the MBH and binary 
stars bound to it; and we find that the probability of ejecting HVSs is substantially enhanced by 
multiple encounters between the MBH and binary stars at a distance substantially larger than their 
initial tidal breakup radii. Assuming that the HVS progenitors are originated from the two thin disks, 
our simulations show that the distributions of the HVS inclination relative to the disk planes can be 
well reproduced by either the mechanism of tidal breakup of binary stars or the mechanism of ejecting 
HVSs by a hypothetical binary black hole (BBH) in the GC. However, an isotropical origination of 
HVS progenitors is inconsistent with the observed inclination distribution. Assuming that the detected 
HVSs were ejected out by tidal breakup of binary stars, its velocity distribution can be reproduced 
if their progenitors diffuse onto low angular momentum orbits slowly and most of the progenitors 
were broken up at relatively large distances because of multiple encounters. Assuming that the HVSs 
were ejected out by a BBH within the allowed parameter space in the GC, our simulations produce 
relatively flatter spectra at the high-velocity end compared to the observed ones; however, the BBH 
mechanism cannot be statistically ruled out, yet. Future deep surveys of HVSs and better statistics 
of the HVS spatial and velocity distributions should enable to distinguish the ejection mechanisms of 
HVSs and shed new light on the dynamical environment surrounding the central MBH. 
Subject headings: Black hole physics — Galaxy: center — Galaxy: halo — Galaxy: kinematics and 
dynamics — Galaxy: structure 



1. INTRODUCTION 

Hypervelocity stars (HVSs) discovered recently are 
mainly B type stars with mass Z-AMq and ra- 
dial velocities up to 700 km s~^ escaping away 
from the Gala ctic halo (Brown et al."2005'; 'Hirsch et alj 
20051 lEdelmann et all [2005; Brown et al. 2007, 2009aF 
Bromlev et al.ll2009| ). These HVSs are believed to be 



ejected out of the Galactic center (GC) by dynamical 
interactions of (bina ry) stars with the c entral massive 
black hole(s) (MBH; iBrown et all l2009al lb[) . There are 
mainly three mechanisms that can produce such high- 
velocity stars: (1) the tidal breakup of binary stars 
in the immediate vicinity of the central MBH (here- 
after TBK mechanism; [H ills 198i jYu fc Tremaindl200l 
iBromlev et "aIll2006HSari et al.ll2010l ): (2) three-body in- 
teractions of single stars with a hypothetical low-mass 
ratio binary blac k holes (BBHs) in the GC (hereafter, 
BBII mechanism; Yu fc Tremainel2003HGua.landris et al 



2005 



20061: 



LevinI 12006 : iBaumgardt et al. 20061 : ISesana et al 
Lockmann fc Baumgardtl 12008 ): and (3) interac 



tions of single stars with a cluster of stellar mass black 
holes in the yicinity of the central M BH (hereafter, SBHs 
mechanism; lO'Learv fc Loebl I2008D . ^ 

^ Alternative models, other than the GC origin, are a lso pro- 
posed to explain the HVSs by [H eber et al. (2008), Aba di et all 
pOOgD . rriTtukov & Fedoroval l[200a V and .Wang fc Han. ((2009l V 



To determine which mechanism should be responsi- 
ble for the detected HVSs is an important issue in 
understanding the stellar dynamics around the cen- 
tral MBH. Some observational properties of HVSs, 
such as the binarity, rotational velocity, metallicity, 
ejection rate/observed frequency, spatial and veloc- 
ity distribution, etc., are proposed to be helpful in 
distinguishing these mechanisms (e.g., ILu et al.l 



Qistmguishmg these mechanisms (e.g.. ILu et al.l 
HansenI 120071: ISesana et al.l l2007t IPrzvbilla et al.l 



20081 : 



Lopez-Morales fc BonanosH2008l: fPeretsll2009al lb[). How- 
ever, one uncertain factor in these mechanisms is where 
the initial progenitors come from. Different origins of 
the HVS progenitors may also result in different distri- 
butions of HVS properties, especially their spatial dis- 
tribution and the velocity distribution. The main aim 
of this paper is to address this problem through Monte 
Carlo numerical simulations under the constraints from 
current observations. 

The spatial distribution of the detected HVSs is com- 
patible with both the TBK and the BBH mechanism 
if the progenitors of these HVSs are originated from 
disk-like structures in the vicinity of the central MBH 
(e.g., the clockwise rotating young stellar (CWS) disk 
in the GC). The preliminary results that demonstrate 
the disk(s) ori gination o f the H VS progenitors have been 
summarized in lLu et al] ()2010t ) . We shall further demon- 
strate in this paper that an isotropic distribution of HVS 
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progenitors is inconsistent with the spatial distribution 
of the detected HVSs if they are originated from the GC. 
However, the disk(s) origination of the detected HVSs is 
consistent with the distribution of the inchnation angles 
(relative to the disk(s)) of the detected HVSs wh i ch fur - 
ther strengthen the conclusions made in lLu et "ahl ()2010( ). 

The velocity distribution of HVSs is related to not only 
the product ion mechanism but also the origin of their 
progenitors. iSesana et al.l (|2007[ ) have studied the veloc- 
ity distribution of HVSs. They found that the velocity 
distribution of HVSs produced by the TBK mechanism 
for unbound injecting stellar binaries seems to be consis- 
tent with the then detected HVSs though with limited 
statistics, while the HVS velocity distribution produced 
by the BBH mechanism appears to be too flat in com- 
parison with the observations. Their results suggest that 
the HVS velocity distribution may be useful in distin- 
guishing the ejection mechanisms. In this paper, we shall 
further investigate the effects on the velocity distribution 
of HVSs due to different origins of the HVS progenitors, 
e.g., those (binary) stars initially unbound to the MBH 
but later injected into the immediate vicinity of the MBH 
due to some unknown perturbations, and those (binary) 
stars initially bound to the MBH but later evolved onto 
highly eccentric orbits and migrated into the immediate 
vicinity of the MBH. 

The paper is organized as follows. In Section [31 we 
first summarize the observational results on the spatial 
and velocity distribution of the detected HVSs. In Sec- 
tion [3l we explore the detailed dynamics of interactions 
between binary stars on bound orbits and a central MBH. 
The consequences of these interactions are different from 
that between the unbound binary stars on parabolic (or 
hyperbolic) orbits and the MBH intensively invest i gated 
in the hter aturc fe.g.. lHillj[l988t iBromlev et al.l 120061: 
ISesana et a l. 2007). The reason is that the stellar bi- 
nary may experience multiple close encounters with the 
MBH in the former case, while it only experiences a sin- 
gle close encounter in the latter case. Assuming realistic 
distributions of the properties of the initial stellar bina- 
ries, we then simulate both the spatial distribution and 
the velocity distribution of HVSs produced by the TBK 
mechanism and compare the numerical results with the 
observations in Section IH^ In Section (U we also ex- 
plore the interactions between single stars on bound or- 
bits with a hypothesized BBH in the GC. These single 
stars are assumed to be injected into the immediate vicin- 
ity of the BBH from disk-like stellar structures (e.g., the 
CWS disk) surro unding the BBH, wh ich is different from 
that adopted in ISesana et al.l (|2007t) . With reasonable 
but simple assumptions on the parameters of the hypo- 
thetical BBH, the spatial and velocity distributions of 
the ejected HVSs are obtained. Comparison between the 
simulation results and the observations are also discussed 
in Section [SJ The conclusions are given in Section [51 

2. SPATIAL AND VELOCITY DISTRIBUTIONS OF THE 
DETECTED HVSS 

Surveys of HVSs have detected 16 HVSs unbound to 
the Galactic halo, 8 bound HVSs, and 4 HVS candidates 
(|Brown et al.l 120051: iHirsch et al1l2005t lEdelmann et al.l 

^ The cases of interactions between unbound bi naries and the 
centra l MBH have been extensively studied by IBromlev et al.l 
ll2006i) . 



120051: iBrown et ahl l2007l [lOOli)- We summarize their 
spatial and velocity distributions in this section. 

2.1. Spatial distribution 

The spatial distributi on of the HVSs detected so far 
is prob a bly anis otropic (|Abadi et al.ll2009l : iBrown et aH 
I2009bf) . iLu et al . (2010) use great circles to fit the spa- 
tial distribution of the detected HVSs projected on the 
sky of an observer located at the GC, and they find that 
the distribution can be best fitted by two great circles. 
Their results suggest that the spatial distribution of the 
detected HVSs is consist ent with being located on the 
planes of two thin disks (|Lu et al.l 120101 ): (1) eleven of 
the unbound HVSs (plus four bound ones and two can- 
didates; totally 17 objects) are spatially associated to a 
thin disk plane with an orientation almost the same as 
that of CWS disk located within half a parsec from the 
central M BH (see iLevin fc Beloborodov) l2003t iLu et al.l 
[2009: Paumard et al.ll2006l : iBartko et al.lf2009l 120101 ): (2) 
four of the unbound HVSs (plus three bound ones and 
two candidates; totally 9 objects) are spatially associated 
to a thin disk plane with an orientation similar to that 
of the northern arm of the minispiral (or also the outer 
warped part of the CWS disk) in the GC. The normals 
of the best-fit disk planes for these two HVS populations 
arc (?,&)=(311°, -14°) and (176° -53° ) in Galactic co- 
ordinates, respectively (iLu et al.l |2010[ ). Hereafter, we 
refer to those detected HVSs associated with the above 
two best-fit planes as the first population and the second 
population of HVSs, respectively. We denote the inclina- 
tion angle of each HVS to its corresponding best- fit plane 
by Q and describe the spatial distribution of the HVSs 
by a normalized cumulative distribution function of their 
inclination angles P(> 8) (hereafter, 0CDF), which rep- 
resents the number fraction of the HVSs with inclination 
angles higher than 8. The observational 8CDFs for both 
populations of the HVSs are shown in Figure [1] and will 
be compared with the distributions obtained from nu- 
merical models in Sections |4] and [5l For each population, 
we shall compare the P{> 8) of all the HVSs (includ- 
ing unbound HVSs, bound HVSs, and HVS candidates) 
instead of only unbound ones, because (1) for the first 
HVS population, our Kolmogorov— Smirnov (K-S) test 
finds a likelihood of 0.94 that the unbound HVSs and all 
the HVSs are drawn from the same 8CDF; (2) for the 
second population, the number of the unbound HVSs is 
only 4 and the error due to Poisson noise in the 8CDF 
is substantial, therefore we do not show their 8CDF in 
Figure [T] (and uCDF in Figure [2] below, either). 

The gravitational potential of the Galaxy is not exactly 
spherical, and its non-spherical component may deflect 
the radial trajector ies of HVSs after t hey were ejected 
from the GC Ce.g.. lYu fc Madaul 12001 . Given the dis- 
tance and the velocity span (30 kpc < R < 130 kpc 
and 690 km s'^ < v < 980 km s-\ see Section 1^ of 
the detected HVSs, however, the deviation due to the 
flattening of the Galactic disk and the triaxiality of the 
Galactic halo is at most severa l degre es (typically 2°) as 
demonstrated in lYu fc Madaul (j2007t ). Compared with 
the standard deviation of 8 (^ 7° — 8°) shown in Fig- 
ure [1] that small deviation can be ignored. 

2.2. Velocity distribution 
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The detected HVSs have been decelerated in the Galac- 
tic gravitational potential after they were ejected from 
the GC. To easily compare with results obtained from 
theoretical models and see the dependence of HVS ve- 
locity distributions on different ejection mechanisms, we 
first remove the velocity deceleration (caused by the 
Galactic bulge, the disk, and the halo; see below) from 
their observed values with correction for the proper mo- 
tion and then obtain their velocities at infinity by assum- 
ing that they move from the GC to the infinity only in 
the gravitational potential of the central MBH. We de- 
note these HVS velocities at infinity by i;^. and their 
CDF P(> v"^) (hereafter, wCDF) represents the number 
fraction of HVSs which have velocities at infinity higher 
than v°f. 

Generally, the Galactic gravitational potential can be 
described by four components, $ = $BH+<l'buigo+^disk+ 
^haio, i-C-, the contribution from the central MBH, the 
Galactic bulge, the disk, and the halo. We have $bh = 
—GM,/r, where the mass of the centra l MBH in the GC 
is adopted to be M. = 4 x 10^ Mp, (e.g.. lGhez et al.ll2008l : 
iGillessen et al.|[2009t) . In this pap er, we adop t the m odel 
for the last three components from lXue et alj (|2008[ ). i.e.. 



'&bulgc 



GMbuigo 
r + r-buigo ' 



(1) 



with Mbuigc = 1.5 X 10 A/q and the core radius rbuigc = 
0.6 kpc, 

GMdi.,(l - e--/") 

5'disk = , [^) 

r 

with Mdisk = 5 X 10^"A/q and the scale length & = 4 kpc, 
and 



$halo = - 



^^%^ln(l + ^), (3) 

with = ^liere pc is the cos- 

mic critical density, Avir = 200, il^ is the cosmic frac- 
tion of matter, the virial radius rvir = 267 kpc, and 
the concentration c = 12. To see the velocity deceler- 
ation of a star caused by the bulge, disk and halo poten- 
tials above, the minimum velocity required for the star 
starting from the potential center (in the absence of the 
MBH) is - 645 km s~i (or 670 km s"\ 700 km s'^) so 
that they can move to a distance of 50 kpc (or 100 kpc, 
200 kpc) further away. With the Galactic potential 
model above, we correct the velocity deceleration and 
obtain the z;CDFs for both HVS populations, which are 
shown in Figure [21 

In the estimation of above, we adopted the spher- 
ical approximation of the Galactic disk potential given 
by Equation (2), since all the detected HVSs are mov- 
ing at Galactocentric distances much larger than the 
disk sca le length b. If using the Miyamoto— Nagai po- 
tential (jMivamoto fc Nagail 119751 ) to describe the flat- 
tened potential of the Galactic disk, i.e., $(i?, z) = 

-GMdisk/(i?' + (adisk + (z' + 6Lk)'^')')'^' and Afdisk = 
5 X lO^Af0, the scale length and height of the disk are 
fldisk = 5 kpc, fodisk = 0.4 kpc, respectively, we find that 
the difference in the estimated v°? due to the disk flat- 
toning is negligible (^ 10 km s"-'^). 

Adopting a different Galactic potential model may 
result in a different wCDF. For example, choosing 



th e Galactic pot e ntial as that given by Equation (8) 
in iKenvon et al.l ()2008l ). the resulted ranges from 
850 km s~^ to 1200 km s~^, which is substantially 
higher than that shown in Figure [2] (from 690 km s~^ 
to 980 km s~^). However, the slope of the vCDF in this 
velocity range is only slightly f latter than that estimated 
from the pote ntial given by (IXue et al.l l2008t also see 
discussions in ISesana et al.ll2007() . Therefore, our con- 
clusions made in Sections |4] and [5] on comparisons of the 
simulated uCDF (mainly the shape) with the estimated 
uCDF here are not sensitive to the choice of the Galactic 
potential (see further discussions in Section 143]) . 

Note here that the Galactocentric distances of the de- 
tected HVSs range from 25 kpc to 130 kpc, which corre- 
spond to their travel time from the GC ranging from a 
few tens to 260 Myr. 

3. TIDAL BREAKUP OF BINARY STARS BOUND TO THE 

MBH 

The progenitors of the detected HVSs are probably 
originated from disk-like stellar struct ures in the GC 
(e.g., the CWS disk) as demonstrated bv lLu et al.l()2010D . 
These progenitors are likely to be bound to the central 
MBH. In order to investigate whether the statistical dis- 
tribution of the detected HVSs shown in Section [5] can 
be reproduced by the TBK mechanism, here we first ex- 
plore the detailed physics of the tidal breakup of binary 
stars that initially bound to the central MBH as it is 
still elusive in the literature. Consider that the binary 
star is moving on an eccentric orbit bound to the MBH 
and injecting/migrating into the vicinity of the MBH. 
For convenience, we denote the injecting/migrating stel- 
lar binary by the "stellar" binary with semimajor axis Ob, 
eccentricity eb, and mass m = mi + 7712 (where mi and 
m2 are the masses of its two components, respectively). 
The center of mass of the stellar binary and the MBH 
also compose a binary system, which we denote by the 
"outer" binary with semimajor axis Oout and eccentricity 

^out ■ 

3.1. Three-body experiments 

We perform three-body experiments of dynamical 
interactions between stellar binaries and the cen- 
tral MBH. We use the code D0RPI5 based on 
the explicit fifth(fourth)-ordcr Runge— Kutta method 
(jDormand & Prince 1980; Hair er et all |1993[) to inte- 
grate the set of the equations that control the three-body 
interactions between the stellar binary and the MBH. 
We set the tolerance of the fractional energy error to be 
10~^ for a single orbit integration. We have also checked 
that statistically our results are not affected by setting 
a smaller minimum tolerance of the energy error. The 
initial conditions set for the stellar and the outer bina- 
ries include (1) the masses of the two components of the 
stellar binary (mi and m2, m2 < mi) and its mass ratio 
q — m2/mi; (2) the initial semimajor axis and eccentric- 
ity of the stellar binary (ab.i, eb,i); (3) the initial semima- 
jor axis and eccentricity of the outer binary (aout,i, eout,i); 
(4) the orbital orientations of both the stellar and outer 
binaries; and (5) the initial orbital phases of both the 
stellar and outer binaries. 

For convenience, we also define a dimensionless pen- 
etration parameter D to describe how close a stel- 
lar binary can approach the MBH relative to its tidal 
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(degree) © (degree) 

Fig. 1. — Observational inclination distributions of HVSs relative to their best-fit planes (0CDF, P(> 0)). The left panel is for the 
first population of HVSs which are spatially associated writh a disk plane writh an orientation similar to the CWS disk, and the right panel 
for the second population of the HVSs which are spatially associated with a disk plane with an orientation similar to the northern arm 
(Narm) of the minispiral or the warped outer part of the CWS disk. The solid histograms represent the distributions of all the detected 
HVSs, for each corresponding population, including the unbound HVSs, the bound HVSs and the HVS candidates. The distribution of the 
unbound HVSs is also shown for the first population by a dashed line. See Section l2.1l 
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Fig. 2. — Observational cumulative distribution of the ejection v eloc ity at infinity for the two populations of the detected HVSs (i)CDF, 
P(> v^)). Other legends arc the same as in Figurc[T] See Section 12.21 



breakup radius rtb ( see a similar definition in lHillslll98 
iBromlev et aIll2CI06h . i.e., 



D 

where r. 



100^ 

rth 



Ob.i 



3M. 



10'5(toi + TO2) 



-1/3 



(4) 



flout, i(l — Gout.i) is the initial pericenter dis- 
tance and rtb = [3il/,/(mi -|- m2)]^/^ab,i. The penetra- 
tion parameter D = 100 corresponds to the tidal breakup 
radius of a stellar binary, and the range of D of interest 
in this paper is D < 250 where tidal breakup of stellar 
binaries is possible. 

3.1.1. Initial settings 

For the numerical experiments in Sections 13.21 and 13.31 
we simply set the mass ratio of the stellar binary as q = 1 , 
and mi = TO2 = 3Mq; but for those in Section |4] we 
adopt more realistic distributions of the mass ratio and 
the masses of the two components of the stellar binaries 
according to current observations on binary stars. 

We set the initial semimajor axis ab,i of the stellar 
binary in the range from 0.05 AU to 2 AU.^ Pertur- 
bations on a binary star with ab.i < 0.05 AU from the 

^ The timescale for those stellar binaries to be disrupted 



tidal field of the MBH easily lead to the merger of its 
two components; while tidal breakup of binary stars with 
ab.i > 2 AU usually leads to ejections of stars with veloc- 
ities substantially smaller than the hypervelocities inter- 
ested in this paper. For demonstration only, we simply 
choose aba = 0.1 AU and eba = (or 0.1, 0.3, and 0.6 
alternatively) in the three-body experiments presented 
in Sections (3 21 and [33t but in Section SI we adopt a dis- 
tribution of ab,i based on the constraint obtained from 
current observations and eb,i = 0. 

We set the semi major axis of the outer binary a s 
Oout.i ~0.04-0.5pc ()Lu et al.ll200llGillessen et al.ll2009t ). 
as the stellar binaries are possibly originated from the 
CWS disk within half a parsec from the central MBH. 
The pericenter distance of the binary is set to be close 
to the tidal radius of the stellar binary. In Sections 13.21 
and 13. 31 we set the penetration parameter D typically to 
be ~20— 250 and choose aout.i = 0.2pc for demonstration 

by their interac tions with backgroun d stars is roughly 10^ yr 
II Yu fc Tremaine| [2003: Hop manI 120091) . If these binaries are ro- 
tating around the MBH with semimajor axis of < 1 pc as adopted 

in this paper, the total time for ~ 1000 orbits is roug hly < 10'' yr. 
Therefore, the evolution of stellar binaries due to their dynamical 
interactions with background stars can be safely ignored. 



Spatial and velocity distributions of HVSs 



5 



only. 

The orientation of the orbit of the stellar binary, rela- 
tive to that of the outer binary, (f> g [0, tt], is assumed to 
be uniformly distributed in cos 4>, if not specified. 

The initial orbital phases of both the stellar and outer 
binaries are randomly chosen in all the following calcu- 
lations. 

3.1.2. Approximations 

As the stellar binaries are initially set on orbits weakly 
bound to the central MBH, binaries with large penetra- 
tion parameters (e.g., D > 150) may revolve around the 
central MBH for many (e.g., 10—100 or even 1000) or- 
bits. In our calculations below, the period of the stellar 
binary is usually much smaller than that of the outer 
binary. Therefore, most of the calculation time may be 
spent on integrating the stellar binary orbit when it is 
faraway from the MBH. However, when the binary star 
is faraway from the MBH, the tidal torque from the MBH 
on the binary star is rather weak, so we can approximate 
the motion of the stellar binary around the MBH into two 
independent two-body problems: one is for the stellar bi- 
nary, and the other is for the outer binary on an elliptical 
orbit, both of which can be done analytically. We adopt 
the above two-body approximation when the tidal force 
from the MBH on the stellar binary is less than 10~^ of 
the gravitational force between its two components (we 
have checked the two-body approximation by setting a 
lower threshold, e.g., 10~^ or 10~*, and found no signifi- 
cant difference in our results). With this approximation, 
our calculation time is sped up substantially when Oout is 
large (> 0.2pc) and its accuracy can still be maintained. 

A star may be tidally disrupted if its close passage to 
the MBH is < r*i^ = R^{M,/m^y^^ , where m* and i?, 
are the mass and radius of the star, respectively. The ra- 
dius of a star can be roughly given by i?* oc RQ{m^,/MQ) 
for m* < Mq and i?* cx Rc^jm^, I Mp,)^'^^ o therwise, 
where Rq is the solar radius ()Torres et al.ll20ldl and ref- 
erences therein). Due to the tidal disruption, part of 
the disrupted remnants may be swallowed by the cen- 
tral MBH and part may be ejected out. In our calcu- 
lations, we terminate the evolution of the correspond- 
ing system once any component of the stellar binary ap- 
proaches within a radius r*;^ from the MBH. 

The tidal torque from the MBH may change the semi- 
major axis and the eccentricity of the stellar binary dur- 
ing its close passage to the MBH (see further discussions 
in Section [3.2. ip . In the three-body experiment, we as- 
sume that the two components of the binary star merge 
into a single star once the distance between the two com- 
ponents becomes smaller than the sum of their radii, and 
the evolution of the corresponding system is terminated 
once a merger event occurs. 

The consequences on the dynamical interactions of 
stellar binaries with the MBH can be different depend- 
ing on whether the stellar binaries are initially on bound 
or unbound orbits. If the binary star is initially bound 
to the MBH, it may revolve around the central MBH 
for many orbits before its disruption or merger; while 
the binary initially unbound to the MBH passes by the 
MBH only once even if it is not broken up during its close 
passage (see further discussions in Sections 13.21 and l3.3p . 
Multiple times of interactions of bound binary stars with 



a cent ral MBH have been discussed by lAntonini et al.l 
(|2010l ) , in which they focus on the post-Newtonian effects 
on stellar orbits in the gravitational field of the MBH 
and mergers of the two components of the binary. In 
Section [321 we analyze the consequence of the first close 
encounter of the stellar binary with the central MBH, 
and the changes in the orbital elements of the stellar bi- 
nary (if the binary star survives). Then we consider the 
cumulative effects of multiple encounters in Section 13.31 
We focus on discussing how the spatial and velocity dis- 
tributions of HVSs are connected with their originations, 
ejecting mechanisms, and dynamical environments in the 
GC. 

In Section 13. 2[ IC* three-body experiments are per- 
formed for each of the four set of initial conditions, i.e., 
eb,i = 0, 0.1, 0.3, and 0.6, respectively. If a steUar bi- 
nary is not tidally broken up during its first close pas- 
sage to the MBH, the three-body experiment ends up 
when the stellar binary reaches the apoapsis after its 
first close passage to the MBH. In Section [331 10"* three- 
body experiments arc performed for each set of initial 
conditions and those experiments that do not lead to 
disruption/merger/separation within 500 revolutions are 
excluded. 

3.2. Consequences of the first close passage of a stellar 
binary to the MBH 

Considering that binary stars are injected toward the 
MBH on weakly bound orbits with eccentricity close to 1 
and D < 250, the changes in their orbits, as consequences 
of the tidal effect from the MBH, after their first close 
passages, can be characterized by the following five cat- 
egories. 

1. The two components of the stellar binary remain 
bound to each other, but its orbital parameters (ob, 
Cb) are changed by the tidal torque from the MBH 
during the close encounter. The outer binary is 
still bound and the stellar binary will encounter 
with the MBH again in its next close passage. The 
probability of our three-body experiments below 
resulted in this category is denoted by Pbound- 

2. The binary star is tidally broken up, with one com- 
ponent of the binary being ejected out as a high- 
velocity star, and the other one being trapped onto 
an orbit that is tightly bound to the MBH com- 
pared to its parent stellar binary. The probability 
resulted in this category is denoted by Pcj • 

3. The binary is tidally broken up into two single 
stars, and the two stars are not bound to each other 
but both bound to the central MBH. The proba- 
bility resulted in this category is denoted by Pgp. 

4. The two components of the binary merge into one 
single star. The probability resulted in this cate- 
gory is denoted by Pmrg- 

5. One (or both) component (s) of the binary is (are) 
tidally disrupted and partly swallowed by the cen- 
tral MBH. The probability resulted in this category 

is denoted by Ps^a,llow(= 1 - Pbound - fej - Psp - 
Pmrg)- 
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3.2.1. Consequence 1 

If a stellar binary can still maintain its integrity after 
its close passage to the central MBH, the eccentricity and 
the semi-major axis of the binary are generally changed 
due to the tidal torque from the MBH. 

• Changes in eccentricities. Figure [3] shows the sim- 
ulation results of the changes in eccentricities of 
the stellar binaries after their first close passages 
to the MBH. When the penetration parameter is 
small (e.g., D < 150), the eccentricity of the stel- 
lar binaries may be excited (or de-excited) to var- 
ious values in the range from to 1. The large 
range of the excited (or de-excited) eccentricities 
is mainly due to various phases of the stellar bi- 
naries at the close passages to the MBH and var- 
ious orientations of the stellar binaries relative to 
the orbit of the outer binary. For encounters with 
D > 200, excitation (or de-excitation) on eccentric- 
ities is relatively small, especially for those binaries 
with eb,i = 0. Multiple subsequent encounters of 
the stellar binary with the MBH should be common 
at Z? > 150, where both Pej and Pmrg are small and 
Pbound is large (see Figured]). In these cases, the 
binary eccentricity (and semimajor axis) may then 
be further changed due to the cumulative effect of 
multiple subsequent encounters (see Section [3. 3p . 

• Changes in semimajor axes. As shown in Figure 21 
the semimajor axis of a binary star can be excited 
to a larger value but can also be de-excited to a 
smaller value. The changes of the semimajor axes 
depend on the orbital phases of the stellar binary 
during its close encounter. If the relative posi- 
tional vector of the two components of the stellar 
binary is aligned with the radial vector from the 
central MBH, the stellar binary is more likely to 
be stretched and its semimajor axis is changed to 
a larger value after the encounter; while the binary 
can also be shrunk and its semimajor axis may be 
changed to a smaller value if the relative positional 
vector of its two components is close to perpen- 
dicular to the radial vector. The time for the bi- 
nary passing by the periapsis is comparable to, if 
not substantially smaller than, the orbital period 
of the stellar binary when the closest approach is 
roughly the tidal radius of the binary. Therefore, 
whether the binary is stretched or shrunk is de- 
termined by the orbital phase and orientation of 
the stellar binary during its close encounter with 
the central MBH. As the stellar binary is stretched 
in most of its orbital phases, the semimajor axis 
of the binary is more likely to be excited rather 
than de-excited after its first close passage to the 
MBH (see Figure |4]). We also find that the relative 
change of the semimajor axis of the stellar binary 
is independent of the initial semimajor axis of the 
binary, which is verified by our calculations with 
adopting different ab,i. 

In addition, we note here that the orbit of the outer 
binary star almost remains the same after the first 
close encounter if the stellar binary can maintain 
its integrity. Our calculations above show that 



the semimajor axis of outer binaries may change 
by Saout ^ 0.02aout and the distance to the peri- 
center of the outer binary may only change by 
Srp < O.OOlrp, which suggests that the location of 
the subsequent close encounters is almost the same 
as that of the first encounter if no other perturba- 
tion takes place on the orbit of the binary star. 

3.2.2. Consequences 2 and 4-' ejection and merging 
probabilities 

Figure [5] shows our calculation results on the ejection 
probability Poj and the merging probability Pmrg- As 
seen from Figure [5l Poj increases with decreasing pen- 
etration parameters D because the tidal force from the 
MBH increases with decreasing D. And Pej also increases 
with increasing initial eccentricities of the stellar bina- 
ries. The reason is that the distance between the two 
components of a binary for given D and ab, and hence 
the tidal torque on them, is larger if the binary has a 
larger eccentricity, and therefore the binary is easier to 
be broken up. For the same reasons, Pmrg shown in Fig- 
ure [5] decreases with increasing D at the large-D end, 
but decreases with decreasing D at the low-D end due 
to the significant increase of Poj . A stellar binary with a 
larger initial eccentricity is also easier to be excited to an 
eccentricity high enough for its two components to merge 
and thus Pmrg increases with increasing eb,i at large-Z? 
end. 

Figure [6] shows the ejection velocity at infinity (v^) of 
the ejected component as a function of the penetration 
parameter (D) of stellar binaries with different initial 
eccentricities (eb.i = 0, 0.1, 0.3, and 0.6, respectively). 
We find that the rms values of the ejection velocities can 
be well fitted by the following formula: 



(KT) 



2\l/2_. 



''cj,0 



+ Gba 
mi -I- 7712 



1-r^y 

260 / 



6Mp 



1/3 



M. 



Qb,i 
0.1 AU 

1/6 



-1/2 



4 X lO^M, 



(5) 



for 20 < £> < 260, where Doff = 13/^1 + e^s is the ef- 
fective penetration parameter for a stellar binary with 
non-zero eccentricity, and v^q = 2560 km s^^ is the 
normalization of the ejection velocity when ChA = and 
ab,i = 0.1 AU. For stellar binaries with zero initial eccen- 
tricities and any D in the range of 20—150, this fitting 
formu la gives a value simila r to that given by Equation 
(1) in lBromlev et all (|2006[) . 

If flout,! is set to a smaller value, the ejection veloc- 
ity may be changed by — ^ ^(u^)^ — GM, /aout,i be- 
cause the ejected star has to first overcome its initial 
bounding energy. If Oout,! = O.Olpc, for example, the 
ejection probability Poj is significantly smaller than that 
shown in Figure [S] and the ejection velocity is also 
substantially smaller because a significant fraction of the 
stellar binaries end up in consequence 3 instead of con- 
sequence 2. 

3.3. Multiple encounters 

3.3.1. The ejection/merging probabilities and the number of 

revolutions 



Spatial and velocity distributions of HVSs 



7 



b,1 



b,1 



0.9 



0.0 



(a) 



(b) 




50 



100 



150 



200 



50 



100 



150 



200 



D 



D 



Fig. 3. — Eccentricity of the binary stars after their first close encounters with the central MBH (e^ i) vs. the penetration parameter 
D, if the binary star survives. The semimajor axes and the masses of the binaries are all set to be aj, ; = 0.1 AU and mi = = 3Mq, 
respectively. The semimajor axes of the outer binaries are all set to be flout,! = 0.2pc. Panels (a)-(d) show the ei, i for those stellar binaries 
with initial eccentricities of e^, ; = 0, 0.1, 0.3, and 0.6, respectively. As seen from this figure, the eccentricity of a binary may change to 
various values between and 1 after the first close encounter if Z3 < 150. The total numbers of survived binary stars in the 10* three-body 
experiments here are 7256, 7160, 6709, and 4782 in panels (a), (b), (c), and (d), respectively. 

The tidal effect from the MBH may be not large 
enough to break up a stellar binary with a large D (e.g., 
D > 150) during the first close passage of the binary. If 
the binary star is initially bound to the MBH, multiple 
times of close passages of the binary to the MBH should 
be common, and thus the probabilities of ejection and 
stellar coalescence are substantially enhanced, especially 
at large D, as demonstrated in Figure [71 which can be 
understood as follows. 



probability is given by 



Qd,0 



1 - Phou 



id,0 



(6) 



The increase of the number of close passages to the 
MBH. Assuming that the ejection probability after 
each close passage of the binary is roughly the same 
as denoted by Poj.o and the probability of the bi- 
nary maintaining its integrity is roughly Pbound,Oj 
then after A'orb revolutions the cumulative ejection 



If iVoib ^ 00, then Pej = Pej,o/(l - Pbound,o) > 

Poj,o- When D is large (> 150), Pej can increase 
substantially with increasing TVorb as Pbound.o is 
close to 1. 

The increase of the semimajor axis of the stellar 
binary. The semimajor axis of the stellar binary 
may be excited to be larger than its initial value 
after the close passages (see Figure (HI), and corre- 
spondingly the penetration parameter D cx Tp/ob 
decreases as rp is almost not changed. A smaller 
D results in a larger Poj. As demonstrated in Fig- 
ure m the semimajor axes of many stellar binaries 



8 



Zhang, Lu, & Yu 



b,1 b,i 




a /a 

b,1 b,i 



1 r-. 



(a) 



(c) 




(b) 




(d) 



•••.•-.••\iv5a~/^ 




50 



100 



150 



200 



50 



100 



150 



200 



D 



D 



Fig. 4. — Relative changes in the semimajor axes of the binary stars after their first close passages to the central MBH as a function of 
the penetration parameter D. The initial conditions for the stellar binaries are the same as that in Figure [3] The initial eccentricities of 
the stellar binaries are i = (panel (a)), 0.1 (panel (b)), 0.3 (panel (c)), and 0.6 (panel (d)), respectively. Generally, the change in a^, is 
large if D is small. The total number of survived binary stars shown in each panel is the same as that in Figure |3] 



before their breakup and last close encounters with 
the MBH are substantially larger than their ini- 
tial values. As the tidal force is approximately a 
monotonically increasing function of Ob, the semi- 
major axes of many binaries tend to rapidly (expo- 
nentially) increase after many subsequent close en- 
counters with the MBH, while the semimajor axes 
of some other binaries decrease only slightly, which 
leads to the asymmetric distribution of ab,f/ab,i, 
i.e., the relative change in ah for many binaries 
with ah excited is substantially larger than that 
for a relatively small number of binaries with Ob 
de-excited, as shown in Figure HI 

• The increase of the eccentricity of the stellar bi- 
nary. Binary stars with large penetration parame- 
ters {D > 150) may be hard to break up initially. 
However, their eccentricities are relatively easy to 
be excited to large values due to cumulative tidal 
effects from the central MBH during their multiple 



encounters with the MBH. The increase in the ec- 
centricity thus leads to a larger ejection or merging 
probability as mentioned in Section [321 We note 
here that the Kozai mechanism may cause periodi- 
cal oscillations of the internal eccentricity of stellar 
binaries during their multiple encounters with the 
MBH, which may be important in enhancing the 
merger of the two components of those binaries. A 
detailed study of thi s resonance can be found in 
lAntonini et"all (j201(]| ). 

For those stellar binaries injected into the vicinity of 
the MBH with D -200-250,' - 90% of them are dis- 
rupted or merged within 500 revolutions. For those stel- 
lar binaries that are broken up within 500 revolutions 
with one component being ejected out. Figure [HI shows 
the cumulative probability distribution of the number 
of revolutions around the MBH (A'orb) that they experi- 
enced. The majority of the stars are broken up within 10 
revolutions when D is small, while most of them can be 
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Fig. 5. — Ejection and merging probabilities (P^j and Pmrg) 
of the stellar binaries after their first close passages to the MBII 
as a function of the penetration parameter D (top panel and bot- 
tom panel, respectively). The initial conditions set for the stellar 
binaries are the same as those in Figure |3] The solid, dashed, 
dotted, and dot-dashed curves represent the results for those stel- 
lar binaries writh initial eccentricities e^, ; = 0, 0.1, 0.3, and 0.6, 
respectively. 

broken up only after many revolutions when D is large. 
If 250 > D > 200, for example, iVorb ranges from 4 to 
500 (see Figure [9]), and about 40% of those stellar bi- 
naries are broken up after 100 encounters. We should 
note that the ejection velocities of the resulted HVSs at 
large D are substantially smaller than those of the stellar 
binaries broken up at small D. 

3.3.2. The effects of relative orientations 

The tidal effects on a stellar binary during its close pas- 
sages to the central MBH are also affected by the initial 
orientation of the stellar binary orbit with respect to the 
outer binary orbit. Given the penetration parameter D, 
the prograde rotating stellar binary is easier to be tidally 
broken up and can receive a larger kick in velocity dur- 
ing the close passages than a retrograde rotating stellar 
binary. Compared to prograde rotating stellar binaries, 
retrograde rotating binaries can revolve around the cen- 
tral MBH for many more orbits before they are broken 
up. And the ejection velocities of the stars from the in- 
teraction of retrograde binaries with the MBH are also 
much smaller than that for the cases of prograde bina- 
ries. In Table [U we list some numerical results to show 
the effect of different relative orientations of the stellar 
binary with respect to that of the outer binary. (For 
demonstration, only the results from those experiments 
with the stellar binaries initially having ab,i =0.1 AU, 
eb,i = 0, and D = 100 are shown in Table [TJ) 



TABLE 1 

Effects of the relative orientation of the stellar binary 

ORBIT TO the outer BINARY ORBIT 



Orientation 




-Pmrg 






CT b 




Prograde 

Retrograde 

Uniform 


0.88 
0.66 
0.86 


0.12 
0.29 
0.13 


0.000 
0.038 
0.004 


2.39 
1.17 
1.80 


+0.13 
-0.15 
+0.53 

+0.27 
-0.37 


1.1 

16 
3.9 



Note. — In the first column, "prograde" (or "retrograde") denotes 
that all the stellar binaries are rotating in the same direction as (or the 
anti-direction of) the outer binary, while "uniform" denotes that the 
orientations of the stellar binaries with respect to that of the outer binary 
(p are uniformly distributed in cos 4> for (f> £ [0, tt] . This table shows the 
results from those experiments with the stellar binaries initially having 
ab,i =0.1 AU, eb,i = 0, and D = 100. 

^ Mean ejection velocity at infinity in units of 10'' km s~^. 
^ Errors in the rms of the ejection velocities at infinity. 

Mean number of close encounters experienced by the stellar binaries 
before their breakup and the ejec tion of one of its two components (i.e., 
the consequence 2 in Section l3.2|l . 



TABLE 2 

Different models for the TBK mechanism 



Model 


/(aft,,)" 


firp) 


lout, 


(Pc) 


D?^i( km s-i)!" 


UB-1 


Opik 


Constant 






250 


UB-2 


log-N 


Constant 






250 


LP-1 


Opik 


Constant 


0.04- 


-0.5 




LP-2 


log-N 


Constant 


0.04- 


-0.5 




RW-1 


Opik 


Random walk 


0.04- 


-0.5 




RW-2 


log-N 


Random walk 


0.04- 


-0.5 





The "log-N" denotes those mode ls adopting a log-normal dis - 
tribution of /(ab.i) as that given bv lDuguennov &: Mavo3 II1991I ). 
and the "Opik" denotes those models adopting the Opik law for 
/{«b,i). 

^ For the "UB" model, those stellar binaries initially have veloci- 
ties of 250 km s~^ at infinity. 

3.3.3. Defl,ection angles 

Here, we briefly summarize the results about the de- 
flection angles between the direction of an HVS moving 
away from the Galactic halo and the initial injecting di- 
rection of its progcnitorial binary. The ejecting direc- 
tion of an HVS is almost anti-parallel to the injecting 
direction of its progenitor and the off-set 5Q is < 15°, 
as demonstrated by iLu et al] (|2010l . see Figure 1). In 
another words, HVSs can well memorize the injecting 
directions of their progenitors. We also find here that 
the major component of 5Q is in the orbital plane of 
the outer binary and the other component of the off-set, 
which is perpendicular to the orbital plane of the outer 
binary and is relatively insignificant ^ 1°— 2°. The main 
reason for this result is that the stretching on the stellar 
binary by the MBH tidal force is most effective when the 
orientation of the stellar binary is similar to that of the 
outer binary orbit. Considering the gravitational poten- 
tial due to the stellar cusp surrounding the MB H (which 
currently seems not to b e clearly determined; (jPo et al.l 
120091: iBartko et ani2010( )). the orbit of a stehar binary 
revolving around the MBH may precess significantly if 
the binary is broken up only after more than hundreds 
of revolutions. Due to this precession, the ejecting direc- 
tion of a resulted HVS may deviate significantly from the 
injecting direction of its progenitor and the main compo- 
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Fig. 6. — Ejection velocity at infinity (fj?) for tfie ejected components of those stellar binaries that are tidally broken up after their 
first close passages to the central MBH. The initial conditions for the stellar binaries are the same as that in Figure [S] Panels (a), (b), (c), 
and (d) show the results for those stellar binaries with initial eccentricities ej, ; = 0, 0.1, 0.3, and 0.6, respectively. The filled circles are 
the square roots of the rms values of (f^)^ in different bins of penetration parameters D, and the error bars indicate the square roots of 
the standard deviations of {v'^)'^ . The solid curves show the best fits to the rms values. The total numbers of the experiments shown are 
2261, 2339, 2687, and 3543 in panels (a), (b), (c), and (d), respectively. 



ncnt of the deviation is in the orbital plane of the outer 
binary. If projeeting these two directions onto the plane 
that is perpendicular to the orbital plane of the outer bi- 
nary, however, the difference between the two projected 
directions is still insignificant. This suggests that the 
spatial distribution of the ejecting HVSs can still main- 
tain on the plane that their progenitors originated even 
if the precession due to the stellar cusp is significant. 
Therefore, the spatial distribution of HVSs should re- 
flect the geometrical structure of the parent populations 
of their progenitorial binaries. If the progenitors of HVSs 
are originated from a disk-like structure, particularly, the 
resulted HVSs projecting to the infinity on the sky should 
locate close to a great ci rcle with the sa me orientation as 
the disk (see Figure 5 in lLu et "aLll2010l) . We will demon- 
strate this further in Section |4l 



4. INCLINATION AND VELOCITY DISTRIBUTIONS OF 
HVSS FROM THE TBK MECHANISM 

4.1. Initial settings 

To investigate the spatial and velocity distributions of 
HVSs, we now choose more realistic distributions of the 
parameters for both the stellar and outer binaries (i.e., 
initial conditions) involving in our numerical simulations 
as follows. 

For the stellar binary, the initial conditions arc set as 
follows under the assumption that they are originated 
from stellar disk(s) like the CWS one, though many im- 
portant issues related to the disk origination of HVSs 
(such as warping and relaxation of the disk(s) and many 
youn g stars ob served in the GC but not in the CWS 
disk: iLu et al] 120091: iBartko et al.ll2010l ) need to be fur- 
ther investigated (|Lu et al .Il2010[) . il) The distribution of 
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Fig. 7. — Probabilities of different consequences of binary stars 
after N^^h times close passages to the central MBH as a function 
of the initial penetration parameter D. The solid, the dotted, and 
the dashed curv es represent P^j, Pmrg, and Pboundj respectively 
(sec Section l3.2l l. For each type of the curves, we show the results 
for N^rb = li 20, and 500, respectively (solid and dotted curves: 
from bottom to top; dashed curves: from top to bottom). 




D 

Fig. 8. — Ratio of the scmimajor axis of the stellar binary to 
its initial value {a^ f/a^.i) vs. the penetration parameter D, where 
Ob f is the scmimajor axis of the stellar binary when it is at the 
apocenter of the last outer binary orbit just before its breakup. 
For those stellar binaries with D > 150, most of them have ab,f 
substantially larger than a^^i. In the figure, the total number of 
the stellar binaries that are broken up is 6007, of which 2116 are 
broken up after their first close passages to the MBH and thus have 
Qb.f/ab.i = 1- 
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Fig. 9. — Cumulative probability distributions of the number of 
close passages that the stellar binaries have experienced before their 
breakup with ejecting one of their two components. Only those 
stellar binaries that are broken up within 500 revolutions around 
the MBH are taken into account. The solid, dashed, dotted, and 
dot-dashed lines represent the distributions with different ranges 
of penetration parameters 25 < D < 100, 100 < D < 150, 150 < 
D < 200, and 200 < D < 250, respectively. 



their semimajor axes ab,i f ollows either the Opik law, i.e., 
P{ah)dah oc dah/ah fe.g.. iKobulnickv fc Frveill2()07D. o r 
a log-normal distribution ( Duauennov fc Mavoil Il99i[) . 

(2) The mass distribution of the primary stars mi fol- 
lows the Miller— Scalo initial mass function, i.e., /(mi) cx 
mj"" and a ~ 2.7 (e.g.. lKroupall2002l ).^ For massive bi- 
nary stars, the distribution of the secondary star or the 
mass ratio q = 7712 /mi can be described by two popu- 
lations: (a) a twin population, i.e., about 40% binary 
stars have q ^ 1 and (b) the rest b inaries, which follow 
a dis t ribution of f(q) ~ constan t (jKobulnickv fc Frverl 
[2007t iKiminki et atl 120081 [2009h . (3) The initial ec- 
centricity of the stellar binary is assum ed to be 0, as 
adopt e d in some prev i ous w orks (e.g., iBromlev et alJ 
(|2006( 1: lAntonini et"all ()2010f )). (4) The orientation of 
the stellar binary orbital plane is chosen to be uniformly 
distributed in coscj) for G [0, tt]. 

For the outer binary, the initial conditions are set as 
follows. (1) The orientation of the outer binary is set 
to satisfy a Gaussian distribution around the central 
planes of the stellar disks that initially hosted the stcl- 
lar binaries with a standard deviat ion of ~ 7°— 13°, (c.f., 
iLu et al.|[2009l : iBartko et ai]|2009t) . The orientations of 
the host disks are assumed to be the same as that of 
the two best-fit plan es, i.e., (^ 6)=(3 11°,-14°) and (176°, 
-53°), respectively (|Lu et al.ll2010t) . (2) The stellar bi- 
naries are initially set on orbits with semimajor axes of 
flout,! ~0.04— 0.5pc, and this range is adopted according 
to the observational extents of the CWS disk. We note 
here that the extents of the structure associated with 
the second plane are not clear. Nevertheless, we assume 
that it is in the same range as the CWS disk. Adopting a 
slightly larger range (say, from 0.04pc to Ipc) does not af- 
fect the results presented in this section and next section. 

(3) The distribution of flout,! of these binaries is assumed 
to follow the same surface density distrib ution as that of 
stars in the CWS disk /(flout,i) Ooutj (|Lu et al.ll2009l : 
IBartko et ani2009f) . (4) The distribution of the initial pe- 
riapsis of the stellar binaries should depend on detailed 
mechanisms leading to the injection. It is not clear how 
bound binaries (in the outer disk region) are delivered to 
the immediate vicinity of the central MBH. Nevertheless, 
the mechanism responsible for the small orbital angular 
momentum of the stellar binary may fall into the two 
extreme categories discussed below. 

• Large perturbations on the orbital angular mo- 
menta of binaries initially rotating around the 
MBH with insignificant eccentricities: in this sce- 
nario, changes of the orbital specific angular mo- 
menta of outer binaries are large compared with 
the specific angular momentum required to ap- 
proach the tidal radius rtb (denoted by Jmax — 
y/2GM^rth) , which leads to an even distribution of 
the distance to the periapsis of the outer binaries, 
i.e., f{rp) oc constant. Therefore, many binaries 
can penetrate into the vicinity of the MBH with 
small D, and some of them are broken up during 
the first close passage but others may be broken 
up or merged after many revolutions. Hereafter 

We have tested that adopting a top-heavy initial mass func- 
tion as that in lBartko et al.l 1120101 ) does not significantly affect the 
resulted spatial and velocity distributions. 
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the model adopting f(rp) oc constant for initially 
bound binaries with the above settings is denoted 
by the "LP" model. The distribution of rp above is 
similar to that adopted for ini tially unbound stella r 
binaries which is discussed in ISesana et al.l (|2007| ) . 
For comparison, we also introduce a model (de- 
noted by "UB" below) in which the stellar binaries 
are init ially unbound to th e MBH, similar to that 
used in ISesana et al.l (|2007D . 

• The perturbations on the orbital specific angular 
momenta of binary stars arc small compared with 
"/max, and f{rp) is significantly different from that 
for the "LP" model. In this scenario, the eccentric- 
ity of an outer binary (or rp ) may increase ( or de - 
crease) slowly. For example, iMadigan et al.l (|2009f ) 
demonstrated that a secular instability of eccentric 
stellar disks rotating around the MBH can gradu- 
ally excite some (binary) stars to extremely eccen- 
tric orbits (e — 7- 1). The time period for this sec- 
ular evolution can be hundred to thousand times 
of the orbital period of the outer binary before the 
breakup or merging of the stellar binary. There- 
fore, the probability for the stellar binary to be 
broken up is substantially even at a large D due 
to the cumulative tidal effect from the MBH be- 
fore the binary can possibly migrate into the region 
with a substantially small D (as also illustrated in 
Figure [6l). Note that the ejection velocity of the 
component escaping away from the MBH is sub- 
stantially small if the binary is broken up at a large 
D rather than at a small D. For these cases, the 
velocity distribution of the resulted HVSs may be 
substantially steeper than that from the LP model 
simply due to the suppression of the number of stel- 
lar binaries with small D. The detailed distribution 
of /(r-p) should depend on the detailed perturba- 
tions. Below we adopt a toy model for this scenario 
to obtain the spatial and velocity distributions of 
HVSs and compare them with that obtained from 
the "LP" model. 

The orbit of the outer binary may be scattered by 
weak encounters of the binary star with other stars 
(in or out of the stellar disk), and its angular mo- 
mentum I J| may also be changed after each revo- 
lution. We approximate this as a one-dimensional 
diffusion process of | J|, and the standard deviation 
of the change of the angular momentum after each 
revolution y/ ((<5| J|)^) = constant | J|. This ran- 
dom walk in \J\ corresponds to a change in rp (i.e., 

Svp = ±S,\/rp/rp^o) after each revolution. The Tp.o 
should be set by the penetration parameter (e.g., 
D = 250) where the tidal effect begins to be impor- 
tant for a stellar binary with the maximum semi- 
major axis considered in this paper (obj = 2 AU). 
The periapsis of some stellar binaries may then dif- 
fuse inward to the vicinity of the MBH slowly. Dur- 
ing this diffusion process, these stellar binaries may 
be broken up at a large D after many revolutions 
and few stellar binaries with small Ob.i can move 
into the very stellar region with small D. We de- 
note this model as the "RW" model. For illustra- 
tion, we set the model parameter ^ = 0.4 AU and 
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Fig. 10. — Distribution of simulated HVSs in the —R plane {R 
is the Galactocentric distance of an HVS). The top and the bottom 
panels show the results obtained from the "RWl" model and the 
"LPl" model, respectively. The filled circles and stars represent 
the detected HVSs that arc spatially associated with the best-fit 
plane with an orientation almost the same as that of the CWS 
disk and the best-fit plane with an orientation similar to that of the 
northern arm of the minispiral in the GC (or the outer warped part 
of the CWS disk), respectively. The dots represent the simulated 
HVSs. The apparent lower boundary is due to the selection effects 
of HVS candidates in observations (i.e., the cutoff at the radial 
velocity = 275 km s~^). The vertical dashed lines are the lower 
and upper boundary of the distance of the detected HVSs from the 
GC. The total number of the simulated HVSs shown in the top 
and the bottom panel are 606 and 449, respectively. This figure 
shows that the ejection of HVSs with velocities > 1000 km s^^ is 
significantly suppressed in the "RWl" model. 

rp^o = (i:'/100)rtb ^ 630 AU with Ob = 2 AU. A 
model with a much larger ^ will be reduced to the 
"LP" model above. 

Tabled] lists the above models and the settings of a few 
related parameters. According to the models (i.e., the 
"UB" model, the "LP" model, and the "RW" model), we 
use Monte Carlo simulations to obtain both the 9CDF 
and the uCDF. For each model, the total number of the 
three-body experiments is 10*. We only record those 
cases in which the masses of ejected stars are in the mass 
range (3Mq, 4Mo) of the detected HVSs if not specified, 
and then calculate both the distribution of the inclination 
angles with respect to the central planes of their parent 
population and the velocity distributions. 

To compare with the observations summarized in Sec- 
tion [21 the selection effects must be carefully considered. 
As shown in Figure [101 niost HVSs arc detected in Galac- 
tocentric distances from 25 kpc to 130 kpc which is partly 
due to the detection limit and partly due to the limit 
in the survey area. In the MMT survey of HVSs by 
I Brown et al.l ()2009al) , the HVS candidates are selected by 
a cutoff in the radial velocity, i.e., v^i > 275 km s~^. To 
account for these selection effects, we adopt the Galac- 
tic potential model listed in Section [2] and simulate the 
radial distribution of HVSs under the assumption of a 
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Fig. 11. — Cumulative distributions of HVS inclination angles 
(©CDF) relative to the two best-fit planes of the detected HVSs. 
The histograms in the top and bottom panels represent the ©CDF 
of the first and second populations of the detected unbound HVSs, 
respectively. The solid and dashed curves represent the numeri- 
cal results from the LP-1 model by adopting the thickness of the 
two disks which the HVSs are originated from to be 12° and 13°, 
respectively (also see Figure I12I I. The dashed curves show the 
©CDFs if the stellar binaries are isotropically distributed rather 
than originated from two disks defined by the two best-fit planes 
(see details in Section 14.211 . 

constant HVS ejection rate for each of the above models, 
which appears to be com patible with current observa- 
tions (jBrown et al.ll2007f ). We only take those HVSs in 
the surveyed area with the radial distance from 25 kpc 
to 130 kpc and radial velocities v^f > 275 km s~-^ as the 
simulated sample, which arc then used to compare with 
the observations. The effect of limited lifetime of the 
ejected HVSs is also considered in our calculations. The 
selection effects are similarly considered for those models 
in the BBH mechanism in Section [5] below. We note here 
that the 8CDF is not affected much but the wCDF may 
be significantly affected by the selection effects. 

4.2. 0CDF; the inclination angle distribution of the 
ejected stars 

Figure [TT] shows the 9CDFs obtained from the "RW- 
1" model. As seen from this figure, the observational 
8CDF for both HVS populations can be well reproduced 
if the HVSs were originated from two disk-like stellar 
structures with orientations the same as the two best-fit 
planes and thickness of ~ 12° and 13°, respectively (also 
see the top panel of Figure fT2|) . Here, the thickness of 
a stellar disk is defined by the standard deviation of the 
inclination angle of stars in the disk from the disk central 
plane. For these two populations of the detected HVSs, 
our K-S tests find the likelihoods of 0.896 and 0.999 that 
the observational 0CDF are drawn from the same distri- 




FlG. 12. — Aitoff-projection for both the simulated HVSs and 
the detected HVSs. The red solid circles, the carmine stars, and 
the cyan stars represent the detected unbound HVSs, the HVS 
candidates, and the bound population of HVSs, respectively. The 
top and bottom panels show the simulated HVSs produced by the 
TBK mechanism and the BBH mechanism (the LP-1 model), re- 
spectively. In the top panel, the open circles and triangles repre- 
sent the simulated HVSs which are originated from the two best-fit 
disk planes with normals of («,b)=(311°, -14°) and (176°, -53°), 
and the thickness of these two originating disks are 12° and 13° , 
respectively; while in the bottom panel, the thickness of the two 
originating disks are 7° and 10° , respectively. The green solid curve 
shows the boundary of the survey area in the Northern hemisphere 
in each panel. This figure shows that the spatial distribution of the 
detected HVSs can be well reproduced by both the TBK mecha- 
nism and the BBH mechanism if the progenitors of the detected 
HVSs are originated from two thin disks with the same orientations 
as the fitted ones. One difference in the resulted spatial distribu- 
tions of these two mechanisms is that a small number of simulated 
HVSs from the BBH mechanism can have larger inclination angles 
relative to their original planes, which is due to that some progeni- 
tors can interact with the secondary BH at a very small separation 
and receive a large kick. 

bution as that obtained from the numerical simulations 
for the "RW-1" model, respectively. All the other mod- 
els (i.e., "LP-1", "LP-2", "UB-1", "UB-2", and "RW-2") 
can reproduce the observational OCDFs for both HVS 
populations by choosing suitable thickness (typically in 
the range of 7° — 13°) for the two disks. The resulted 
8CDFs are mainly determined by the thickness of the 
disks where the HVS progenitors are originated, and also 
weakly depend on the mechanism adopted that leads to 
the injection of stellar binaries into the immediate vicin- 
ity of the MBH. 

Alternatively, assuming that the stellar binaries, i.e., 
the HVS progenitors, were isotropically distributed but 
with other initial settings the same as those in the above 
models, the ejected HVSs should also be isotropically 
distributed in the survey area. We do similar numerical 
simulations and obtain the 8CDFs for those isotropically 
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distributed HVSs, relative to the two best-fit thin disk 
planes of the detected HVSs, as shown by the dashed 
curves in Figure [TTJ Here, a simulated HVS is assigned 
to one of the two populations if it is closer to the best-fit 
plane of that population than to the other plane. Us- 
ing the K-S test, we find the likelihoods of 1 x 10""' and 

0. 016 that the eCDFs of the simulated HVSs are the 
same as those from observations for the first and sec- 
ond population of the detect HVSs, respectively. There- 
fore, we conclude that the detected HVSs are highly un- 
likely to be produced from the tidal breakup of isotropi- 
cally distributed progenitorial binary st ars. This furthe r 
strengthens the conclusion obtained in iLu et all (|2010( ) , 

1. e., the detected HVSs are probably originated from two 
thin disks with orientations similar to the CWS disk and 
the northern arm of the minispiral (or the warped outer 
part of the CWS disk) in the GC, respectively. 

4.3. vCDF: the velocity distribution of the ejected stars 

Figure [13] shows the wCDFs for both the simulated 
HVSs (obtained from different models) and the observa- 
tions. Our numerical simulations show that the vCDF is 
almost independent of the thickness of the disk(s) where 
the HVS progenitors are originated, but it does depend 
on how close the stellar binaries can approach the MBH 
and on the initial distribution of the semimajor axes of 
the stellar binaries. Different models produce quite dif- 
ferent wCDFs. In the "LP" models and the "UB" models, 
for example, more than 30% of the resulted HVSs (with 
> 700 km s"-') have velocities larger than the max- 
imum velocity of the detected HVSs (~ 1000 km s~^), 
while the "RW" models can produce a steep wCDF which 
is quite similar to the observational ones. The models 
with log-normal distributions of Ob.i produce less HVSs 
at the high- velocity end because of the fraction of stellar 
binaries with small Obj (i.e., < 0.3 AU) is substantially 
smaller compared with that in those models with the 
Opik law. For the first HVS population, our K-S tests 
find 2.7 X 10-6 (7.9 x lO'^) and 2.5 x 10"^ (1.3 x lO^^) 
likelihoods that the uCDFs obtained from the "LP-1" 
("UB-1") model and the "LP-2" ("UB-2") model are 
drawn from the same distribution as the observational 
ones, respectively, which suggests that the first HVS pop- 
ulation is unlikely to be produced by either of the "LP" 
model and the "UB" model. For the second HVS popu- 
lation, the K-S likelihood is 0.01 (0.05) and 0.06 (0.07) 
for the "LP-1" ("UB-1") model and the "LP-2" ("UB- 
2" ) model, respectively. These numbers suggest that the 
second population is not likely to be produced by the 
"LP" or "UB" models with an initial Ob.i distribution 
of 1/ab.i but it may not be inconsistent with the "LP" 
(or "UB") models with a log-normal distribution of ab,i 
(though with limited statistics). However, the wCDFs 
resulted from the "RW" models appear to be consistent 
with the observations as the K-S tests find the likeli- 
hoods of 0.52 (0.13) and 0.43 (0.52) that the observa- 
tional uCDFs of the first (second) HVS population are 
the same as that obtained from the "RW-1" model and 
the "RW-2" model, respectively. 

Adopting a different form of the Galactic potential 
may affect the estimation of the wCDF for the detected 
HVSs in Section 2 and the selection effects discussed in 
Section 4.1. For example, if adopting a simple Galac- 
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Fig. 13. — Cumulative distribution of the ejection velocity at 
infinity (dCDF) obtained from different models for the TBK mech- 
anism. The top panel and the bottom panel represent the first and 
second populations of the detected HVSs, respectively. The thick 
(thin) dot-dashed, dashed, and solid curves represent the vCDF 
obtained from the "LP-1" ("LP-2") model, the "UB-1" ("UB- 
2") model, and the "RW-1" ("RW-2") model, respectively. Both 
the "RW-1" and "RW-2" models can reproduce the observational 
vCDF, while the other models generate a significant number of 
HVSs with large that are not detected by current HVS survey. 

tic potential mo d el as that described by Equation (8) 
in iKenvon et al] (|2008D . the of the detected HVSs 

ranges from 850 km s"^ to 1200 km s"'^ and the slope 
of the uCDF is slightly flatter than that estimated in 
Section 2. The simulated tjCDFs from the models of 
"LP-l", "LP-2", "UB-1", and "UB-2" are not likely to 
be consistent with the vCDF of the detected HVSs, while 
both the "RW-1" model and the "RW-2" model can pro- 
duce a uCDF similar to that estimated for the detected 
HVSs according to the new Galactic potential. Our main 
conclusion that the TBK mechanism can reproduce the 
detected wCDF made in this section is not affected by 
th e choice of the Gala ctic potential (also see d iscussion 
in iSesana et all ()2007l ) and lKenvon et al.l ()2008f )). 

To close this section, we note here that the fraction of 
stellar binaries resulting in ejection of HVSs with prop- 
erties similar to the detected ones is around ^ 3%— 6% 
in those models adopted above (see Table [2]). Current 
observations imply that the total n umber of HVSs sim - 
ilar to the detected ones is ~ 100 ()Brown et al.ll2007t) . 
Therefore, the number of stellar binaries is required to 
be around a few thousands and the injecting rate is 
~ lO-^yr"^. This rate appears to be roughly consistent 
with that estimated by iMadigan et al.l ()2009t ) if the pro- 
genitors are injected into the immediate vicinity of the 
MBH due to secular instability developed in the stellar 
disk. However, we caution here that the dynamics lead- 
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Fig. 14. — Parameter space for the hypothetical BBH in the 
GC if it is responsible for the ejection of the detected HVSs. The 
short dashed line represents the upper limit of the current semi- 
major axis of the B BH that i s compatible with oth er observations 
l|Gualandris fc Merritti 120091 : lYufc Tremaini[2003 'l. The vertical 
solid line represents the upper limit of the mass ratio of the BBH 
if the secondary BH now exists in the GC. The dot-dashed line is 
the constraint on the initial semimajor axis of the BBH from the 
HVS ejection velocity (see Section |5. It . The dotted line is the con- 
straint on the initial semimajor axis of the BBH from the longest 
travel time of the detected HVSs, which should be smaller than the 
decay timescalc of th e BB H only by gravitational wave radiation. 
See details in Section |5. II 

ing to the injection ot (binary) stars is currently not clear, 
and therefore any argument based on the event rate may 
have substantial limitation in distinguishing the produc- 
tion mechanism of the detected HVSs. 

5. EJECTING HVSS BY A HYPOTHETICAL BBH IN THE 

GC 

The dynamical interactions between a BBH and single 
stars can also eject HVSs (the BBH mechanism). In Sec- 
tion 15.11 we discuss possible constraints on the possible 
parameter space of a BBH if it exists (or existed) in the 
GC and is responsible for the ejection of the detected 
HVSs. According to those constraints, in Sections [5?l]- 
15. 4[ we perform a large number of Monte Carlo simu- 
lations for a few models of the decay of the BBH orbit 
to investigate whether the spatial and velocity distribu- 
tions of the detected HVSs are compatible with the BBH 
mechanism.^ 

5.1. Parameter space for the hypothetical BBH 

The parameter space for the secondary BH, if be- 
ing currently located close t o the central primary 
MBH, has been investig a ted by [Hansen fc Milosavlicvid 
(I2003D. lYu fc Tremainel^(l2003fi. .Gillessen ct al, 120091 ) . 
and " iGualandris fc Merritti (|2009[) . The observations 
so far have not shown evidence in contradict with 
the existence of a BBH with semimajor axis aeBH 
2mpc and r nass ratio v = Mt_2 /M,_-\ < 0.02 (see 
Figure 2 in lYu fc Tremaind (|2003f ) and Figure 13 in 
IGualandris fc Merritti l\26m ). where M, i and M, 2 are 



the masses of the two components of the BBH. If a BBH 
is responsible for ejecting the detected HVSs, simple con- 
straints on the BBH parameter space may be directly 
obtained by the properties of the HVSs as follows (see 
Figure [T4|). 

If multiple intermed iate-mass BHs exist in the GC 



UPortegies Zwart et al][20M ). HVSs can also be ejected from the 
GC by interactions with the BHs; but the stability of the system is 
not clear and an A'^body simulation of ejection of HVSs from such 
a system is beyond the scope of this paper. 



• The travel time of the detected HVSs ranges from 
30 to 260 Myr so that the BBH should exist in 
the GC 260 Myr ago and should not have merged 
into a single BH 30 Myr ago, which put a strin- 
gent constraint on the location of the hypothetical 
BBH in the GC. The timescale for the decay of the 
orbit of a BBH in the gravit ational wave radiation 
dominated stage is given by (|Peterslll964[ ) 

rr "BBH 
J^GW = 



|aBBH| 
5 



'BBH 



64 G3 Af.,iA/.,2(A/..i+M.,2) 

= 3.6 X 10 V 



\ mpc 



r'(e) 



M,, 



4 X IO^Mq 



where ctBBH 
/(e) ^(1 



)^(l-e^) 



da^mildt and 
73 



-7/2 



1 



24 



37 

96^ 



(7) 



(8) 



Combining Equation ([7]) and the constraint 
from the HVS travel time, the hypothetical 
BBH should have a semimajor axis obbh.o ^ 

2.2 (z^ + z^^)^^'' mpc when the detected HVSs be- 
gan to be ejected out from the GC (see the long- 
dashed line in Figure [T4|) . Note that the above 
constraint on cbbh.o is obtained by assuming that 
the initial eccentricity of the BBH is bbbh.o = 
(/(gbbh.o) = 1) and the BBH orbital decay is dom- 
inated by the energy loss due to gravit ational wave 
radiation (for details see lPeterslll964|) . Adopting 
a moderately large initial BBH eccentricity does 
not change the constraint much (e.g., obbh.o is in- 
creased only by a factor of 2.5 if bbbh.o = 0.5). 
However, the constraint can be changed by orders 
of magnitude if the initial eccentricity is close to 1. 

• Before the gravitational wave radiation domi- 
nated stage, the orbital decay of the BBH when 
it is "hard" with obbh.o < G'M,,2/4ct^ ~ 
430i/(cr/100 km s-i)-2mpc (|Quinlanlll99^ should 
be dominated by the interactions of the BBH with 
unbound stars in the Galact ic b ulge or boun d stars 
in th e GC stellar cusp (e.g., |^i2002; Sesan a et al.l 
I2OO8I ). And the interactions between the BBH and 
the stars injected from the stellar disk (e.g., the 
CWS disk) may be insignificant to the BBH or- 
bital decay as the number of those stars injected 
into the vicinity of the BBH is relatively small com- 
pared with that from other bound stars. Although 
HVSs can be ejected out as the high-velocity tail 
when the BBH semimajor axis is large (e.g., from 
0.2pc to a few mpc, see ISesana et al.l l2008l and 
iBaumgardt et al.l I2006D . the ejecting efficiency is 
low and the majority of the ejected high-velocity 
stars should be old stars from the cusp rather than 
young stars originated from the stellar disk. If the 
progenitors of the detected HVSs were originated 
from the stellar disk surrounding the MBH and 
the HVSs are ejected out by the BBH mechanism, 
the BBH (with a given mass ratio) should be hard 
enough to efficiently eject such high- velocity stars. 
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The rms velocity of the ejected HVSs at infinity, if 
their progenitors are on parabolic (or hyperbolic) 
orb its with eccentricities close to 1, can be given 
by (|Yu fc Tremainell2003D 



TABLE 3 

Different models for the BBH mechanism 



740 km s" 



0.01 



/ aBBH 

\ nipc 



-1/2/ M.1 



4 X lO^Mn 



(9) 

Therefore we should have a rough constraint on the 
BBH semimajor axis when the BBH began to eject 
those detected HVSs, i.e., aBBH,o ^ llOi/mpc, as 
the minimum velocity of these HVSs at infinity is 
^ 700 km s~i (the solid line in Figure [T?|) . Consid- 
ering that HVSs can be ejected by the BBH with a 
larger semimajor axis, we reset aBBH.o < 220i/mpc 
which corresponds to an rms ejection velocity of 
^ 500 km s^i and is substantially smaller than the 
value at which the BBH becomes "hard" ([Quinlanl 
[l996h . 

5.2. Models for the orbital evolution of the BBH 

The orbital evolution of a BBH during its hard 
stage is dominated by three-body interactions of low- 
angular momen t um st ars with the BBH. As shown in 
lYu fc Tremaind (I2003D . the decay timescale for the hy- 
pothetical BBH in the GC is roughly 4-8 Gyr during 
the hard stage unless the potential of the Galactic bulge 
is significantly flattened or triaxial (jYul I2002D or there 
are a large number of massive perturbers ( Perets et al.l 
l2007t) which lead to efficient transferring of stars onto 
low-angular momentum orbits. This decay timescale is 
substantially longer than the travel time of the detected 
HVSs (< 260 Myr). According to these estimates and 
the constraints on the parameter space of the BBH, we 
can use the following three simple models to generally 
describe the orbital evolution of the BBH during the 
course of ejecting the detected HVSs. A few parameters 
involved in these models are listed in Table [H 

(a) If the potential of the Galactic bulge is spherical or 

at least is not significantly fiattened/triaxial and 
the massive perturbers in the Galactic bulge is not 
sufficient in quickly filling the loss-cone, then the 
orbital decay of the BBH is slow. In this model, we 
assume that the semimajor axis of the BBH does 
not change during the period of ejecting the de- 
tected HVSs and aBBH,o — 220z^mpc. (Choosing a 
somewhat smaller value does not affect our results.) 

(b) If the BBH is initially at cbbh.o = 

2.2 (i^ 4- 1^2) ^^^mpc and the orbital decay is dom- 
inated by gravitational wave radiation, we have 
the evolution of the semimajor axis of an initially 
circular BBH as obbh = aBBH.o(l ~ t/260Myr)i/"'. 
For a BBH with a moderately large initial eccen- 
tricity eBBH,0: the orbital evolution is similar to 
the circular one. For a BBH with extremely larg e 
initial cbbh.o (e.g., see iMatsubavashi et al.l 120071 ) . 
the decay timescale is too short and thus the BBH 
should not be able to eject HVSs over a time span 
similar to the longest travel time of the detected 
HVSs. 
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^ The origin of the injecting stars is set to be from a disk similar to the 
CWS disk for the first five models listed here, while it is set to be from 
infinity and unbound to the BBH for the last model (the "FIX-UB" 
model). 

(c) If the orbital decay of the BBH during the hard 
stage is much faster than that of the spheri- 
cal case because o f a highly fiattened or triax- 
ial bulge 1^ [2001 or many m assive perturbers 
in the bulge (jPerets et al.l I2007D . the semimajor 
axis of the BBH may quickly evolve from the up- 
per boundary to the lower boundary shown in 
Figure [141 In order to account for the detected 
HVSs, this decay timescale should be comparable 
to the span of the HVS travel time. Otherwise, 
it is either reduced to model (b) or simply incon- 
sistent with observations. In this model, we as- 
sume an exponential evolution for the BBH semi- 
major axis, i.e., obbh = obbh.o cxp(— t/260Myr) ~ 
220i'exp(— t/260Myr)mpc, which corresponds to a 
constant aBBH/|aBBH|- 

With the above models, we perform Monte-Carlo sim- 
ulations of the three-body interactions between a single 
star and the hypothetical BBH with mass ratios v = 0.01 
and 0.003, respectively (see Table |3l). For each model, 
the total number of the three-body experiments is 10'*. 
We set the tolerance of the fractional energy error to be 
10~^ for each three-body experiment. The initial condi- 
tions for the injecting stars are set to be similar to those 
given in Section |4| for the outer binaries. The initial dis- 
tribution of the penetration parameters, now defined as 
'"p/obbh.Oi is set to be /(^p/aBBH.o) constant (simi- 
lar to the "LP" models for the TBK mechanism), and 
the starting point is '''p^o/aBBH.o ^1.1—1.2 because the 
ejecti on efficiency and ve locity drop rapidly for rp > rp.o 
(e.g.. ISesana et al.ll2006l ). The orientation of the BBH 
relative to the normal of the stellar disk(s) is set to a few 
typical values, e.g., 6* = 0, 7r/4, 7r/2, etc., and our calcu- 
lations show the resulted 9CDF and uCDF are slightly 
affected but not much by the choice of this orientation. 
We can also choose the "RW" models for the BBH mech- 
anism (similar to that for the TBK mechanism) , but our 
calculations show this model do not lead to a suppression 
of the ejection of HVSs at the high- velocity end. 

5.3. eCDF 

Figure[l5]shows the GCDFs obtained from the "GW2" 
model listed in Table |3| for the BBH mechanism. As seen 
from the figure, the observational OCDFs can be well re- 
produced by this model if the HVS progenitors are orig- 
inated from the two best-fit disks with thickness of 7° 
and 10°, respectively (also sec the Aitoff-projection of 
the simulated HVSs in the bottom panel of Figure IT^ . 
For these two populations of the detected HVSs, our K-S 
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Fig. 15. — Legends are the same as in Figure [TOl but the sohd 
curves represent the ©CDF obtained from the "GW2" model in the 
BBH mechanism. The sohd curves in the top and bottom panels 
show the 0CDF of the simulated HVSs which arc originated from 
the two disk-like structures with orientations the same as the two 
best-fit planes and thickness of 7° and 10° , respectively. 

tests find the likelihoods of 0.97 and 0.91 that the obser- 
vational 0CDF are drawn from the same distributions 
as that obtained from the numerical simulations for the 
"GW2" model, respectively. All the other models listed 
in Table [3] can produce GCDFs similar to the observa- 
tional ones by choosing suitable thickness of the two disks 
(typically in the range 7° — 10°) where the HVS progeni- 
tors are originated. As seen from Figure [T2l a number of 
HVSs resulted from the BBH mechanism can have sub- 
stantially large inclination angles relative to their original 
planes; while no such HVSs appears in the TBK mecha- 
nism. The reason for this difference is that there always 
exist some HVS progenitors which can interact with the 
secondary BH at a very small separation and receive a 
relatively large kick. These very close interactions lead 
to the large inclination angles of some simulated HVSs. 
This difference may be helpful in distinguishing the HVS 
production mechanism in future. 

The deviation of the simulated HVSs from the central 
plane of the disk(s) (with a given thickness) is larger if 
the mass ratio of the BBH is larger (see the analysis in 
iLu et al]l2010t ). But for the allowed parameter space of 
the hypothetical BBH as shown in Figure [T31 the differ- 
ences in the resulted 9CDFs are quite small for different 
choices of the mass ratio (and the semimajor axis) and 
they are easily compensated by slightly different choices 
of the disk thickness. Similarly, a different choice of the 
orientation of the disk plane relative to the orbital plane 
of the BBH may also result in slightly different 0CDF, 
but which can also be compensated by different choices 
of the thickness of the disk(s). 
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Fig. 16. — Cumulative distribution of the ejection velocity at 
infinity (dCDF) obtained from different models for the BBH mech- 
anism. The histograms in the top and bottom panels represent the 
observational iiCDFs for the first and second populations of the 
detected HVSs, respectively. The solid, long-dashed, short-dotted, 
dot-dashed, dot-dot-dashed, and dotted curves represent the vCDF 
obtained from the "FIX" model, the "GWl" model, the "GW2" 
model, the "GW3" model, the "ExpD" model, and the "FIX-UB" 
model, respectively. 

5.4. vCBF 

FigurcfTOlshows the vCDFs for both the detected HVSs 
and the simulated HVSs (obtained from different evo- 
lutionary models for the BBH mechanism listed in Ta- 
ble [3]). Our calculations show that the uCDF (with 
> 700 km s~^) is almost independent of the choice 
of the thickness of the disk where the HVS progenitors 
are originated. The simulated f CDF is insensitive to the 
choice of the orientation of the disk relative to the BBH 
orbital plane. The simulated vCDF is also insensitive 
to the choice of the mass ratio v (if v ~ 0.01—0.001) 
and the eccentricity of the BBH (also see ISesana et alJ 
I200I . As seen from Figure 1161 different models produce 
quite similar t;CDFs and too many HVSs with ejection 
velocities substantially higher than those of the detected 
HVSs. The primary reason for these similarly flat wCDFs 
is as follows. The HVS progenitors injected into the re- 
gion with penetration parameter rp/aeBH ^ 1 can gain 
some energy during the dynamical interaction with the 
BBH, and the mean energy gain is mainly determined 
by the semimajor axis and the m ass ratio of the BBH 
(|Quinlanlll996l : ISesana et al.ll2006[ ) but insensitive to the 
detailed values of rp/oBBH- With increasing Tp/aBBH to 
be slightly larger th an 1, the energy g ain dramatically 
drops to (also see ISesana et al.l [20061 ). Therefore, the 
slope of the uCDF is largely determined by the scatters of 
the energy gains around their rms value at ^p/aBBH ^ 1, 
and the scatters are mainly due to the different orbital 
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phases of the BBH and insensitive to other BBH pa- 
rameters. Because of these features in the BBH mech- 
anism, even the "RW" model for injecting stars, which 
suppresses the ejection of HVSs at the high-velocity end 
for the TBK mechanism, does not lead to such a sup- 
pression for the BBH mechanism here. 

With the relevant parameters for the models listed in 
Table [3l the simulated wCDFs are all quite flat (except 
for the "FIX-UB" model), compared with the observa- 
tional ones. For the first HVS population, our K-S tests 
find that the likelihoods for the simulated wCDFs to be 
consistent with the observed distribution are 0.01, 0.18, 
0.003, 0.0001, 0.0001, and 0.01 for the "FIX", "FIX- 
UB", "GWl", "GW2", "GW3", and "ExpD" models, 
respectively. Except the high likelihood for the "FIX- 
UB" model, the other small K-S likelihoods above sug- 
gest that the first population of the detected HVSs are 
unlikely drawn from the same distributions as those ob- 
tained from the other models. For the second HVS popu- 
lation, the K-S likelihoods arc 0.03, 0.12, 0.05, 0.02, 0.02, 
and 0.04 for the "FIX" , "FIX-UB", "GWl", "GW2", 
"GW3" , and "ExpD" models, respectively. Those likeli- 
hoods are not small enough so that the wCDF for the sec- 
ond HVS population does not appear to be inconsistent 
with those obtained from the those models, especially for 
the "FIX-UB" model. For the combined sample of both 
the first and the second HVS populations, the KS likeli- 
hoods are 0.03, 0.14, 0.003, 3 x 10-^ 3 x 10-^ and 0.02, 
respectively. 

For the models in which obbh evolves with time (i.e., 
the "GW" and "ExpD" models), our simulations indi- 
cate a weak dependence of the ejection velocity on time, 
as the BBH semimajor axis is smaller at a later time and 
the corresponding ejection velocity is larger (see Equa- 
tion (9)). Co mpared to ou r resul ts, the simulated vCDF 
obtained by ISesana et al.l (|2007( ) is flatter, as in their 
model the BBH decays fast at large separations and more 
stars with higher velocities are ejected at a later time as 
the BBH becomes harder. Due to the limited statistics 
of the detected HVSs and the uncertainties in the BBH 
dynamical evolution model, it appears to be premature 
to use the wCDF to distinguish whether the HVSs are 
ejected due to the TBK mechanism or the BBH mech- 
anism and also constrain detailed dynamical models for 
injecting stars to the MBH(s) vicinity. 

Note that the ejected stars are actually moving in the 
Galactic potential, some of those with low velocities (e.g., 
< 400 km s~^) may return to the vicinity of the BBH 
and receive further kicks, which may change the final 
uCDF. We check this effect by doing additional three- 
body experiments for the ejected stars with low ejection 
velocities. We find that the obtained final wCDFs are 
only slightly flatter than those shown in Figure [16] and 
our conclusions above will not be significantly affected. 

6. CONCLUSIONS AND DISCUSSIONS 

In this paper, we first use three-body experiments to 
study the interactions of the MBH with binary stars 
bound to the MBH. We find that the probability of eject- 
ing HVSs is substantially enhanced due to multiple en- 
counters between the MBH and the stellar binaries in- 
jecting into its vicinity even at a distance substantially 
larger than the tidal breakup radius. Given a penetra- 
tion parameter (e.g., D > 150), the velocities of the 



HVSs ejected after multiple encounters are substantially 
smaller than those ejected by first encounters because of 
the excitations of their semimajor axes and eccentricities. 

Assuming that the progenitors of the detected HVSs 
are originated from stellar disk structures, by using 
Monte Carlo simulations we find the distribution of the 
inclination of HVSs relative to the disk planes can be 
well reproduced by both the mechanism of tidal breakup 
of binary stars and the mechanism of ejecting HVSs by 
a hypothetical BBH in the GC. However, an isotropical 
origination of HVS progenitors is inconsistent with the 
observed inclination distributions. We find that the spa- 
tial distribution of HVSs is primarily determined by the 
geometrical structure(s) that their progenitors originated 
but less sensitive to whether the ejection is due to tidal 
breakup of binary stars in the vicinity of the central MBH 
or dynamical interactions of stars with a BBH. The se re- 
sults strengthen the conclusion in iLu et al] ()2010l) that 
the detected HVSs are probably originated from the two 
disk-like stellar structures in the GC, one of which is 
probably the CWS disk. 

We find that the HVS velocity distribution should en- 
code some dynamical information in the GC environ- 
ment. Assuming that the detected HVSs were ejected 
out from the GC by tidal breakup of binary stars, its ve- 
locity distribution can be reproduced if the HVS progen- 
itors diffuse onto low angular momentum orbits slowly 
and most of the progenitorial binaries were broken up at 
large distances with small ejection velocities. In this sce- 
nario, the HVS velocity distribution not only depends on 
the distribution of the initial semimajor axes of the stel- 
lar binaries but also on how the binary stars diffuse onto 
low angular momentum orbits to reach the vicinity of 
the central MBH. If the progenitors were injected to the 
vicinity of the MBH by large perturbations on their or- 
bital angular momenta, the simulated HVS velocity dis- 
tributions appear to be fiattcr at the high-velocity end, 
which are inconsistent with the observed ones. If the 
HVSs were ejected out by a BBH in the GC, our sim- 
ulations show that the HVS velocity distribution is not 
sensitive to the mass ratio and other properties of the 
BBH but ma.y dep end on how the BBH orbit decays (cf., 
ISesana et al.ll2007l ). We find that it is less likely that the 
observed HVS velocity distribution can be reproduced 
with the allowed parameter space of the BBH, as the 
BBH mechanism produces a relatively fiat spectrum at 
the high-velocity end; however, the BBH mechanism can- 
not be statistically ruled out, yet. Future deep surveys of 
HVSs and better statistics of the HVS spatial and veloc- 
ity distributions should enable to distinguish the ejection 
mechanisms of HVSs and shed new light on the dynam- 
ical environment surrounding the central MBH. 

So far, most of the detected HVSs are B-type stars, and 
their progenitors are associated to some specific stellar 
structures, such as the CWS disk in the GC. However, 
there are many late-type stars in the stellar cusp of the 
GC and injection of these latc-typc (binary) stars into 
the vicinity of central MBH (or BBH) should result in an 
old population of HVSs in the Galactic halo. The spa- 
tial distribution of the old population of HVSs should be 
more isotropic as the parent population of their progen- 
itors is dynamically relaxed. The velocity distribution 
of this old-population HVSs could also be quite differ- 
ent, as their progenitors may have a different dynami- 
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cal environment and are injected into the vicinity of the 
central MBH differently. Searching for the old popula- 
tion of HVSs (iKollmeier fc Gouldl [20071: iKollmeier et"all 
|2009|) should be helpful in revealing the dynamical envi- 
ronment in the GC. 

Note added in proof: The unbound HVS, HE 0437- 
5439, was proposed to be ejected out from the LMC 
(Edelmann et al. 2005); and this object was excluded in 
the original fits to the two disk planes in Lu et al. (2010). 
However, a recent measurement of its proper motion by 
the HST suggests it should be originated from the GC 
(Brown et al. 2010). As seen from Figure 12, the de- 
viation of HE 0437-5439 away from the disk plane(s) is 
relatively large compared with other HVSs, and it is dif- 
ficult to have its progenitor originated from the best fit 
disk(s) under the TBK mechanism. However, the devi- 
ation can be compatible with the disk origin under the 
BBH mechanism in which a small fraction of HVSs with 



large scatters can be produced (see discussion in Section 
5.3). Because of the short lifetime and long traveling 
time of HE 0437-5439, it is suggested that this object 
was originally a binary star when being ejected out from 
the GC and lately merged into a single blue straggler 
(Brown et al. 2010, Perets 2009b). This proposal is also 
compatible with the BBH mechanism in which hyperve- 
locity binary stars can be ejected out as predicted in Lu 
et al. (2007). 
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